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We discuss selected large-scale anomalies in the maps of temperature anisotropies 
in the cosmic microwave background. Specfically, these include alignments of the 
largest modes of CMB anisotropy with one another and with the geometry and direc- 
tion of motion of the Solar System, and the unexpected absence of two-point angular 
corellations especially outside the region of the sky most contaminated by the Galaxy. 
We discuss these findings in relation to expectations from standard inflationary cosmol- 
ogy. This paper is adapted from a talk given by one of us (GDS) at the SEENET-201 1 
meeting in August 201 1 on the Serbian bank of the Danube River. 
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When we look up at the sky with microwave eyes, we see a uniform glow, such 
as would be emitted by a black body of just under 3°K. This is of course the Cosmic 
Microwave Blackbody (CMB) radiation originating from when the expansion of the 
Universe had cooled and diluted the primordial plasma to the point where neutral 
atoms were finally long-lived, and the Universe thus became transparent. This CMB 
was first measured in the 1960s by Penzias and Wilson. Its discovery was key in 
establishing the Big Bang theory of cosmology's primacy over its rival, the Steady 
State theory. In the 1970s, a dipolar deviation from the isotropy of that signal was 
first measured, reflecting predominantly, it is believed, a Doppler-induced signal due 
to the Solar System's motion through the rest frame defined by the CMB. It took until 
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the 1994 for the COsmic Background Explorer (COBE) satellite's Differential Mi- 
crowave Radiometer (DMR) to finally measure and make a full-sky map of intrinsic 
anisotropics in the radiation at the level of a few tens of ^K. These anisotropics repre- 
sent principally fluctuations in the temperature/density/gravitational potential at the 
time of the photons' last scattering, also known as the epoch of recombination. They 
are thus the seeds of cosmic structures that would grow to become galaxies, clusters 
of galaxies, and, eventually more complicated non-linear inhomogeneites such as us. 

While many experiments in the 1990s and 2000s were able to improve on as- 
pects of the COBE DMR measurements - increasing the S/N and the angular reso- 
lution over small patches - it was only with the publication of the first results from 
the Wilkinson Microwave Anisotropy Probe (WMAP) in 2003 that there was an im- 
proved full-sky map of the CMB temperature fluctuations. One might well wonder 
how one can make a full-sky map of the CMB when the Galaxy intervenes and is 
brighter than the CMB over the part of the sky close to the Galactic plane. Indeed, 
in COBE maps released to the public, the Galactic plane is immediately evident. 
Meanwhile publicity maps from WMAP show no Galactic plane. This is because 
with improved S/N and more frequency bands, WMAP is able to produce a map out 
of linear combinations of the individual waveband maps in which the Galactic fore- 
grounds should be, and at least appear to be, subdominant. This is the so-called Inter- 
nal Linear Combinations (ILC) map. Similar maps have been made by independent 
groups of researchers and are quite close to the ILC. Nevertheless, one should not 
be surprised should one discover remnant Galactic contamination in such maps, and 
for most purposes CMB analysis is done on partial skies in which the Galactic Plane 
and other contaminated regions (such as those surrounding bright point sources) have 
been removed. 

Although the WMAP temperature maps most obviously improve our ability to 
infer the small-angular-scale properties of the CMB, they also, especially because of 
their increased signal to noise, significantly increase our knowledge of the large scale 
properties of the CMB. It is on these large scale properties that this paper focuses. 
In the rest of this paper we will investigate the large angle properties of the CMB. 
We will discover that if one uses the full-sky ILC map then one finds very odd corre- 
lations in the map, that correlate unexpectedly to the Solar System. (Unexpectedly, 
because we were led to believe that any residual contamination in the map would be 
connected to the Galactic foregrounds.) By masking out the region of the sky behind 
the Galactic plane, we will then discover that the sky appears to have remarkably 
little correlation between pairs of points separated by more than 60°, at a level that 
is statistically very surprising. Looking into this anomaly more deeply we will find 
that it remains robust through all seven years of published WMAP data, and further- 
more that it is very difficult to explain within the context of the canonical Inflationary 
Lambda Cold Dark Model model of cosmology. 
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1. MULTIPOLE DECOMPOSITION 




Fig. 1 - The Internal Linear Combination Map is a weighted linear combination of the five WMAP 
frequency maps. The weights are computed using criteria which minimize the Galactic foreground 
contribution to the sky signal. The resultant map provides a low-contamination image of the CMB 
anisotropy. Courtesy of the WMAP Science Team. 



When confronted by WMAP's all sky map of the CMB temperature fluctua- 
tions (Figure [j}, the immediate response of a cosmologist is to expand the map in 
spherical harmonics: 



AT 

it 1 



E ^2 a em Y em (9,(p), (1) 



t=2 m= 



the monopole (£ = 0) and dipole {1 = 1) having been subtracted. This expansion is 
so automatic because the inflationary model tells us that the at m are (realizations of) 
Gaussian random variables of zero mean. (Or nearly so. Non-linear effects can in- 
duce small, but in-principle measurable non-Gaussianity.) The a£ m are therefore the 
most convenient physical variables for comparing observations with theory. More- 
over, in the standard Lambda Cold Dark Matter (ACDM) model, the Universe is 
statistically isotropic, so that the expectation values of at m obey the relation 

(a} m ae' m >) = Ce5ui5 mm > . (2) 

This means that, in the standard theory, the only thing worth measuring is Cg. The 
variances, Q, of the underlying Gaussian variables, ae m , are also the expected values 
of the measured angular power spectrum, 

Q = ^t E W 2 - (3) 

m=— I 
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We shall very sloppily use the same symbol for both the Ci appearing in equation 
[2] and the Ce appearing in equation |3j even though the former is a property of the 
underlying statistical distribution of which the sky is a realization, and the latter is a 
property of the actual sky. More pedagogically careful treatments are readily found 
in the literature, but this level of sloppiness is standard. 

The WMAP angular temperature-temperature (TT) power spectrum is shown 
in Figure [2j (To be specific, this is the power spectrum produced by the WMAP 
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Fig. 2 - The TT Angular power spectrum. The points are the 7-year temperature (TT) power spectrum 
from WMAP. The curve is the ACDM model best fit to the 7-year WMAP data. The plotted errors 
include instrument noise, but not the small, correlated contribution due to beam and point source 
subtraction uncertainty. The gray band represents cosmic variance. Figure is from [ 1 ] courtesy of the 

WMAP Science Team. 

Science Team using the first seven years of data.) Fitting cosmological parameters to 
the inflationary ACDM model allows us to infer important properties of the Universe 
within the context of that canonical model. In particular, given our interest in the 
largest scale properties of the Universe, we learn that the Universe has a geometry 
that is indistinguishable from flat. (Meaning that no curvature, either positive or 
negative can be discerned.) This can be seen approximately by the location of the 
first peak in the power spectrum, which is in the bin centered around £ = 91. (The 
first clear detection of the first peak was made by the TOCO experiment |2], and then 
by the Boomerang collaboration [3 ] who were the first to conclude that the Universe 
is the therefore close to flat.) 

Traditionally, the geometry of the Universe was the only ultra-large scale prop- 
erty of the Universe that one needed to measure; which would lead us to ask whether 
there is anything else interesting to learn about the Universe on largest scales. As we 

submitted to Romanian Journal of Physics ISSN: 1221-I46X 



The Oddly Quiet Universe: How the CMB challenges cosmology's standard model 



5 



shall see, the data suggests that there is. We might begin to suspect that this is the 
case by looking at the angular power spectrum and observing that the value of the 
quadrupole C2 is anomalously small - well outside the grey cosmic variance band. 
Just how unlikely this is has been a matter of extensive, but rather uninformative, de- 
bate. Uninformative, because it is really not the smallness of the quadrupole that we 
shall conclude is really strange about the large-angle CMB sky. Nevertheless, it was 
one of the motivating factors behind various investigators explorations of the \cw-l, 
or large angle, properties of the CMB. 

As this is not a review, we shall not attempt to be exhaustive or even compre- 
hensive in our exploration of large angle CMB anomalies. There are reviews of the 
subject which attempt to be so. Two very different viewpoints are offered by the 
WMAP Science Team itself in Q, and this collaboration in (5J. We shall instead 
focus here on two results that we have highlighted which reflect the corner that the 
large angle CMB seems to paint us into. This is just a small part of the big picture, 
and we apologize to our many colleagues whose fine (!) work we do not cite here. 

First, we shall work with the full sky CMB as represented by the ILC. (It 
doesn't seem to much matter which year ILC one uses, or whether one instead uses 
one of the other full sky maps, so we shall just stick with the ILC map which was 
produced using three years of WMAP data.) We shall look only at the two lowest 
interesting monopoles i = 2,3. In Figure [3] we plot the quadrupole i = 2 plus oc- 




Fig. 3 - Quadrupole plus octopole anisotropy of the WMAP sky map in Galactic coordinates, shown 

with the ecliptic plane and the cosmological dipole. Included are the multipole vectors (solid 
diamonds); two for the quadrupole (red diamonds) and three for the octopole (green diamonds). We 
also show the four normals (solid squares) to the planes defined by the vectors that describe the 
quadrupole and octopole temperature anisotropy. Figure from 1 5 1 . 

ISSN: 1221- 146X 



http: // ww w. nipne . ro/rj p 



6 



Glenn D. Starkman et al. 



topole anisotropy of the ILC sky map in Galactic coordinates. (The £th multipole is 
just Ylm=-e a tmXim-) Various other environmental quantities are shown - the plane 
of the ecliptic (the plane of the Solar System) together with the north and south eclip- 
tic poles (the normal directions to that plane), and the cosmological dipole direction 
(and its antipode). 

We have also included the multipole vectors that describe the quadrupole and 
octopole. Multipole vectors are the analogs for i > 1 of the dipole vector. We nor- 
mally think of a pure dipolar real function f\ (6, (ft) in terms of a vector d, f(9, <ft) = 
d ■ f instead of as a sum of spherical harmonics. (Here f is the unit coordinate vec- 
tor, f = (sm#cos</>,sin#sin</>,cos#).) The magnitude of the dipole is d = \d\ and 
its direction is the unit vector d = d/d. The dipole strength d plus the two degrees 
of freedom of d replace the three real coefficients Re(an), Im(an) and aio of the 
spherical harmonic expansion. (Since the function / is real, a±o is real, and ai_i is 
determined by an.) 

Similarly, we can replace the five real degrees of freedom of the quadrupole 
(ci2 m as constrained by reality conditions) with two unit vectors u 2,1 and u 2,2 and a 
single scalar A^ 2 \ (This is because an angular momentum 2 object can be obtained 
from the product of two angular momentum one objects.) These unit vectors are the 
multipole vectors of the quadrupole. Similarly the seven real degrees of freedom of 
the octopole {a^ m ) can be replaced by three unit vectors u 3 ' 1 , -u 3,2 and -u 3,3 and a 
single scalar A^ 3 \ (This multipole vector representation appeared as early as (6).) 

The multipole vectors of the quadrupole appear in Figure [3] as red diamonds, 
those of the octopole as green diamonds. They are plotted in both northern and 
southern hemispheres because they are defined only up to a sign that can be absorbed 
into A^K We also plot the normals to the plane defined by the two quadrupole mul- 
tipole vectors r^ 2,1,2 ) || (u 2 ' 1 x u 2,2 ) as a red square (again in both hemispheres), 
and as green squares the normals to the three planes defined by the three multipole 
vectors of octopole. Note that the normals cluster together on the sky, implying that 
quadropole plane and the three octopole planes are nearly aligned. Moreover, the nor- 
mals are near the ecliptic plane, implying that not only are these four planes aligned 
but the are nearly perpendicular to the ecliptic. Furthermore the normals are near the 
dipole, meaning that the planes are not just aligned and perpendicular to the ecliptic 
but oriented perpendicular to Solar System's motion through the Universe. Finally, 
as one can see from Figure|3] the great circle of the ecliptic plane (black curve), very 
carefully separates the strong extrema to its south from the weaker extrema to its 
north. The precise statistical significance of these correlations, first discussed in (7J 
(although the alignment of the octopole and quadrupole with one another was first 
pointed out in [8]), depends on how one calculates them. 

However one does the statistical analysis, these apparent correlations with the 
Solar System geometry are puzzling. They do not seem to reflect the Galactic con- 
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tamination that we might have expected from residual foreground contamination in 
the ILC map. Indeed there are a number of challenges to explaining these results. For 
one, the observed quadrupole and octopole are aligned (appear as Ynn in a frame with 
the z-axis along the common or average axis of the four planes), and not as a Yto (in 
any frame). This makes it difficult to explain them in terms of some localized effect 
on the sky. Also, the quadrupole is much smaller than the octopole, which means that 
perturbative explanations in terms of small vectors or gradients are challenging. The 
best one can say is that these full-sky solar-system correlations remain unexplained. 
The Planck experiment will hopefully shed new light on these mysteries. 



2. ANGULAR CORRELATION FUNCTION 

We now wish to leave the harmonic domain and look instead at the real space 
two-point correlation function of the CMB temperature map, 

C(0) = r(rai)r(n 2 )| ftl . ft2=cosfl , (4) 

where the overbar indicates an average over all pairs of directions on the sky sepa- 
rated by the angle 9. 

While it was once traditional to calculate the angular two-point correlation 
function, this has mostly fallen out of favor. Partly this is because of the lore that 
C{9) contains exactly the same information as the angular power spectrum Ct, 

C(0) = ^2 2 ^C £ P e (cos9). (5) 

l 

However, this relationship between observed quantities holds only when evaluated 
on the full sky. If we impose a Galaxy cut on the map before evaluating the Ct or 
C{9), as is typically done, then the relationship fails. Also, it holds in the statistical 
ensemble only if the assumption/prediction of statistical isotropy is correct. It may 
be, but it certainly should be tested. Finally, there is a reason why we often look at 
both a function and its Fourier (or other appropriate) transform - features that are 
hidden in one sometimes become much more evident in the other. 

In figure|4]we plot various versions of the two point angular correlation function 
from the WMAP 5-year data. The smooth dotted line with the blue band around it is 
the C{9) that would be obtained from equation [5] using the angular power spectrum 
Ct predicted by the best-fit ACDM model. That blue band is the one-sigma cosmic 
variance band. In other words, if we vary each of the Cg inside the one-sigma cosmic 
variance range around its expected value, then the C{9) obtained will remain entirely 
within the blue region. 

While the smooth blue band is the theoretical expectation, all of the jagged 
curves are obtained from the data in one way or another. Our first observation is that 
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Fig. 4 - The angular two-point angular correlation function from the WMAP 5 year results. C(6) is 

plotted for maps with their monopole, dipole and Doppler quadrupole subtracted. The V-band 
(dashed-dotted-dotted line), W-band (dashed-dashed-dotted line), ILC (KQ75, dashed line) have had 
the KQ75 mask applied. The full-sky ILC result (solid line) is also shown. Also plotted are C(9) from 
the WMAP maximum likelihood Cg (dotted-dashed line), the WMAP pseudo-C^ (dotted line) and the 
best-fit ACDM Cg. The shaded region is the one sigma cosmic variance bound on the standard ACDM 

theory. Figure from (5). 

none of those data curves look like the theory curve. They do not remain inside the 
blue band. However, one cannot really compare curves by eye because the differ- 
ent points on the C(ff) curve are highly correlated; we must devise some statistical 
measure of the their difference. Nevertheless we shall not focus on the difference 
between the data curves and the theory curve. 

What is actually more striking (and more significant) is that all of the C(8) 
curves that are calculated excluding the region inside the Galactic plane remain re- 
markably close to C(8) = over a very wide range of 9, from about 60 to 170 
degrees. We can quantify this in terms of a statistic first suggested by the WMAP 
Science Team [9]: 

S 1/2 = J ^ dcos9[C(6)} 2 (6) 

We have shown that the p- value of Si/ 2 for the five year ILC outside the KQ75 
Galactic cut is a remarkably tiny 0.025%. 

We see in figure[4]that the six data-derived curves divide neatly into two classes. 
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Four that hug the C{9) = axis closely on scales 9 > 60°, and two that do not. All 
four of the zero-huggers are calculated by taking the sky average in ([5]> only over pairs 
of points both of which are outside the Galaxy (as defined by the WMAP Science 
Team's KQ75 Galaxy cut). In other words they involve direct calculations of C{6) 
over the parts of the sky that are to be trusted. Instructively, the zero-huggers include 
all three individual cut-sky (ie. KQ75-masked) waveband maps - Q, V, W. It is no 
surprise that the cut-sky ILC rounds out the group of four, since the ILC is just a 
linear combination of the individual band maps. 

The other two data-derived C(9) curves also agree with one another, but not 
with the zero-huggers (nor with the theory curve or error band). They include the full- 
sky ILC, i.e. the ILC with no Galaxy cut, and the curve generated by substituting the 
WMAP-reported angular power spectrum Ce (derived using a Maximum Likelihood 
Estimator - MLE) into equation ([5]). Considering the full-sky ILC vs. the cut-sky 
ILC, we learn that essentially all of the large-angle (9 > 60°) angular correlation on 
the sky is due to pairs of points at least one of which is inside the Galaxy, ie. inside 
the part of the sky that we don't trust. 

The difference between the MLE curve and the cut-sky curves is less straight- 
forward, since, as far as we can tell, the WMAP-supplied MLE Ce also derive from 
cut sky data. However, the discrepancy may well trace to the same cause as identi- 
fied in flO} which addresses the argument of [ 1 1 1 that one should first reconstruct the 
full sky from the cut sky and then calculate the full-sky two-point conelation func- 
tion. Reference fT0| shows that in practice the reconstruction of Cf is biased due to 
leakage of information from the region obscured by foregrounds to the region used 
for the reconstruction. This leakage comes because of the need to smooth the map 
before reconstructing. In the region oustide but near the cut, he smoothing incorpo- 
rates data from inside the cut. Since the cut is imposed because the data inside it is 
unreliable, one must decide what to substitute for that data. The results then depend 
on the choice of how to fill the cut. Not surprisingly, the problem is largest when the 
particular Ci one is reconstructing is anomalously small (compared to Ce of nearby 
£). Of course, this is precisely the case for the quadrupole. Reference JTOl did not 
extend the analysis to C{9), or more specifically S\/2, but one may reasonably sus- 
pect that the reconstruction is similarly (maybe even particularly) poor at maintaining 
C{9) = since that property will not be maintained by most choices of how to fill 
the cut region. 

The absence of two-point correlation on large scales is a much larger problem 
than the smallness of C2 that first led the community to worry about the low-£ CMB 
sky. There are two ways to have a sky with low C(9 > 60°). One is to have all the Ce 
small (or at least all the Ce up to some sufficiently high £). This is not our observed 
sky. (For one thing, such a sky would almost certainly also have a low C{9 < 60°).) 
Looking at figure [2] we see that most of the Ce are comparable to their theoretical 
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Fig. 5 - The WMAP first year TT Angular power spectrum. The data are plotted with la 
measurement errors. The solid line shows the best-fit ACDM model from [12]. The gray band around 
the model is the la uncertainty due to cosmic variance on the cut sky. For this plot, both the model 
and the error band have been binned with the same boundaries as the data, but they have been plotted 
as a splined curve to guide the eye. On the scale of this plot the unbinned model curve would be 
virtually indistinguishable from the binned curve except in the vicinity of the third peak. Figure is 
from 1 13 ] courtesy of the WMAP Science Team. 



values. But recall that the WMAP science team used an MLE method to infer the 
low-£ Cg. This was not the case in their first-year data release, so it is better to look 
at figure [5] to confirm that only C2 is particularly small. 

The other way to get C(6 > 60°) ~ (but not C(0 < 60°) ~ 0) is for the 
low-£ Cg to have a particular relationship to one another. What relationship? The 
one obtained by expanding the observed C{0) in the Legendre series of equation [5] 
However, these relationships are delicate and imply that the Cg must be correlated 
with each other. Table[T] taken from [14], shows how, given the higher £ Cg observed 
by WMAP, it is necessary to tune the contributions to C(9) from C2, C3, C4 and C5 
against those from Cq and above in order to get Si/ 2 to be as low as it is (1152 (/iK) 4 ). 
By contrast, in the best fit theory it would have been enough to tune just C2 and C3. 

It is extremely difficult to arrange for the Cg to have particular relative values 
in the context of the standard inflationary model, because the ag m are independent 
Gaussian random variables. Thus, even if we were able to adjust the theory so that 
the expected values of the Cg (i.e. of (21 + l) -1 Ylm \ a em\ 2 ) were precisely what 
was needed to get C(6 > 60°) = 0, cosmic variance would perversely obliterate the 
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Table 1. 

Si j2 (in (/iK) 4 ) obtained by minimizing with respect to Cg (for £ in the range 
2 < £ < i?max,tunc and fixing Cg with £ > ^ m ax,tunc)- We show the statistic for the 
best-fit theory and for WMAP, as a function of the cutoff multipole ^max,tune- 

Also 

shown is the 95 per cent confidence region of the minimized S)/2 derived from chain 1 of 
the WMAP MCMC parameter fit. The bottom row gives the measured value of S\ji 



outside the cut - 1152 (/iK) 4 . Table is taken from 1 14 
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carefully adjusted relationship among the Q of the theory would not be preserved by 
the measured values on the sky in a particular realization of the a,£ m . To be precise, 
when we replace the theoretical Ct by the Cg inferred from a Legendre polynomial 
series expansion of the cut-sky ILC C{9), we find that there is less than a 3% chance 
of recovering an SV^ less than or equal to the observed value of SV/ 2 . Moreover, 
most of those 3% achieve the low Su2 by lowering multiple low-i Cg, and so are 
unlike the observed angular power spectrum. The observed sky, at least the part out- 
side the Galaxy cut, seems not to respect the fundamental prediction of the standard 
cosmological model that the ag m are independent random variables. 



3. SUMMARY 

The inflationary ACDM model has many successes. The ability to fit the peaks 
and troughs of the medium and high-^ CMB TT angular power spectrum with just 
a few parameters is remarkable. However, for the lowest few multipoles and the 
largest angular scales, the observations disagree markedly with the predictions of the 
theory. Examining the lowest interesting multipoles (the quadrupole and octopole) 
of the best full sky CMB map, we find that they appear unexpectedly correlated with 
each other. The plane defined by the quadrupole and the three planes defined by the 
octopole are nearly parallel to each other. They are nearly perpendicular to the plane 
of the Solar System (ecliptic). They point essentially at the dipole - the direction of 
our motion through the CMB. Finally, they are oriented (with respect to their shared 
axis) such that the ecliptic carefully separates the strongest extrema in the north from 
the weaker extrema of the south. (Any review of CMB anomalies would include 
multiple other examples, some of which may well be connected to the above.) 
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These deviations from statistical isotropy in our CMB sky have yet to be ex- 
plained, and there are significant challenges to doing so. Because of their multipole 
structure, these deviations are not characteristic of something one would obtain from 
a mis-handling of the Galactic foregrounds, and they are also difficult to obtain from 
a single localized patch of the sky. They are also not easy to obtain by a perturbative 
expansion in small vectors (errors in the dipole, gradients of potentials, etc) because 
the quadrupole is so much smaller than the octopole. 

Arguably one should not trust the part of the sky behind the Galaxy. Examining 
the two-point angular correlation on the sky outside the Galaxy we find that there is a 
marked absence of correlations above 60° angular separation. By one measure (first 
proposed by WMAP), this absence of correlation has a p-value of just 0.025%; in 
other words, it would happen accidentally in the best-fit ACDM model just once in 
4000 realizations. Most troublingly, it suggests that Ce of different £ are not indepen- 
dent. 

This anomaly too has, so far, found no satisfactory explanation. One could 
imagine that non-trivial cosmic topology could induce covariance among Ci (since 
the fundamental domain does not possess a rotational symmetry), however searches 
for non-trivial topology have so far failed to find any [ 15 -flT]]*] These searches al- 



ready extend nearly to the diameter of the last scattering surface. Searches beyond 
the last scattering surface are possible in principle, but so far none have been demon- 
strated to be powerful. 

Other explanations have been proffered that can reduce Si/ 2 (for example [ 19 ]), 
but none that can reduce its expected value to the observed one. The challenge is how 
to induce covariance among the Ce within the context of the inflationary paradigm. 

Future results from the Planck satellite may show these large-angle/low-£ anoma- 
lies to be nothing more than systematic errors in the measurements or analysis of the 
WMAP (and the COBE) team, but unless and until they do these anomalies remain 
the outstanding point of disagreement between the standard cosmological model and 
observations. 
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